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We analyze a number of proposed explanations for spectroscopic anomalies observed in atomic 
hydrogen defects embedded in a solid molecular hydrogen matrix. In particular, we critically evaluate 
the possibility that these anomalies are related to Bose-Einstein condensation (both global and local). 
For each proposed mechanism we discuss which aspects of the experiment can be explained, and 
make predictions for future experiments. 



I. INTRODUCTION AND MOTIVATION 

Quantum solids, where the zero point motion of the 
atoms is greater than roughly 10% of their separation, 
form a fascinating class of materials. A principal question 
with these materials is to what extent they are quantum 
coherent, and under what conditions they can be super- 
solid - supporting dissipationlcss mass flow. Examples of 
quantum solids include 4 Hei~— , solid hydrogen^, Wigncr 
crystals, and atomic hydrogen defects in solid molecu- 
lar hydrogen^—. Here we theoretically study the last 
system, giving a critical evaluation of scenarios of Bosc- 
Einstcin condensation of atomic hydrogen defects. We 
make testable predictions for these scenarios. 

Solid hydrogen, with a Lindemann ratio of 0.18, is the 
only observed molecular quantum crystal. It is a rich sys- 
tem with many rotational order/disorder transitions 1 ^. 
The phenomenology of H2 solids is even more interesting 
when atomic H defects are introduced^—. Here we focus 
on the spectroscopic properties of this system, and how 
they may be related to quantum coherence. 

Our work is motivated by recent experiments in low- 
temperature (T ~ 150mK) solid molecular hydrogen, 
populated with large densities (n ~ 10 18 cm~ 3 ) of atomic 
hydrogen defects. These experiments observe unex- 
plained internal state populations^—. Ahokas et alZ- 
provocatively conjectured that the anomalies may be re- 
lated to Bose-Einstein Condensation (BEC) of the atomic 
defects. Here our goal is to explore and constrain this 
and alternative scenarios. We conclude that global Bose- 
Einstein condensation could not realistically explain the 
experiments, and that local condensation would lead to 
distinct signatures in future experiments. 



II. EXPERIMENTS 

We review Ahokas et al.'sZ- experimental apparatus, 
results, and observed anomalies. 



A. General introduction: physics of atomic 
hydrogen embedded in solid hydrogen 

Hyperfine structure of atomic hydrogen. Fig. [1] 
schematically shows the level structure for a H atom in 



a B = 4.6T magnetic field, similar to that used in the 
experiments of interest^. At these large fields, the levels 
break into two nearly degenerate pairs. Levels within a 
pair are separated by radio frequencies and the pairs are 
separated by microwave frequencies. The electronic spin 
in states a and b is aligned with the magnetic field and 
is anti-aligned in the other states. 
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FIG. 1: Hyperfine level diagram for hydrogen atom in a B = 
4.6T magnetic field, where the mixing angle is 8 = 3 x 10 -3 . 
Vertical axis schematically denotes energy, while horizontal 
axis has no physical meaning. Arrows denote electron (no- 
slash arrow) and nuclear (slashed arrow) spin projections. 

Ahokas et al£ observe that within the solid hydrogen 
matrix the atomic hydrogen's spectra is modified. The a- 
d energy splitting decreases while the b-c energy splitting 
increases by the same amount. 

Spectroscopy of atomic hydrogen. We consider 
three different spectroscopic probes, distinguished by the 
states they connect: nuclear magnetic resonance (NMR), 
electron spin resonance (ESR) , and electron- nuclear dou- 
ble resonance (ENDOR). NMR uses radio waves to cou- 
ple the a and b states, while ESR uses microwaves to 
couple the a and d or b and c states. ENDOR uses a two 
photon transition to couple a and c. All of these probes 
can be used in the linear regime, where they give in- 
formation about the splittings and occupation numbers, 
or in the non-linear regime. As an example of the latter, 
Ahokas et al£> study what happens when they apply high 
RF power, saturating the NMR line. 

Molecular hydrogen. At these temperatures, two 
states of molecular hydrogen are relevant: the "para" 
and "ortho" states. In the "para" or "p-" configuration 
the relative nuclear wavefunction ip is symmetric under 
exchanging nuclei and the nuclear spins form a singlet. In 
the "ortho" or "o-" configuration ip is antisymmetric and 
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the nuclear spins form a triplet. In all cases, the electrons 
are in a symmetric bonding orbital and the electronic 
spins are consequently anti-aligned. 

The true ground state of solid hydrogen is formed of 
P-H2. In practice, however, it is very rare to have a pure 
P-H2 sample. The ortho state is long lived^, requiring 
hours for the concentration to change by 1%. This ortho- 
para conversion can be a source of heating in experiments 
with an energy A/fc^ = 170 K released per molecule. In 
the experiments of Ahokas et alZ the exact quantity of 
0-H2 is unknown, but given the growth technique it is 
likely to be at least ten percent. 

At standard pressure the P-H2 in the solid is highly 
spherical: interaction with neighboring H2 negligibly dis- 
torts the P-H2. Modeling the hydrogen- hydrogen inter- 
actions by their vacuum values quite accurately describes 
quantities such as the speed of soundiS. If sufficient 
ortho- hydrogen is present, orientational ordering transi- 
tions may occur around IK. The models we consider do 
not rely upon any orientational ordering. Depending on 
sample-preparation conditions, either hep or fee crystals 
may be produced^. 

Atomic hydrogen in the solid lattice. The mo- 
tion of atomic hydrogen in molecular hydrogen has been 
widely studied^. Both thermally activated and quan- 
tum tunneling contribute to defect motion, but quan- 
tum tunneling dominates at these low temperatures (the 
two dominant tunneling pathways have energy barriers 
of 4600K and 100K). One motivation for these studies is 
to explore the possibility of Bose-Einstein condensation 
of these defects. This is conceptually related to superso- 
lidity driven by condensation of vacancies. 

Kumada£ argues on the basis of experimental data that 
the exchange reaction H + H2 — > H2 + H is the dominant 
diffusion mechanism at low temperatures. Other tunnel- 
ing pathways are possible, including correlated, collective 
relaxation and "physical" diffusion. 

Ahokas et alZ achieve populations of 50ppm H defects 
in their solid, and argue that H sits at substitutional 
sites. The observed lifetime of these defects was weeks. 
The dominant decay mechanism should be the recombi- 
nation of two hydrogen defects. One therefore expects 
that this rate is determined by the diffusion rate of the 
defects. The long lifetime is therefore inconsistent with 
the diffusion rates predicted by phonon assisted tunnel- 
ing. One possible explanation is the suppression of tun- 
neling by the strain-induced mismatch of energy levels 
between neighboring site o 11 ^ 2 . 

Crystal growth. Ahokas et alZ grow solid hydro- 
gen from a gas of electron-spin polarized metastable hy- 
drogen atoms, which undergo two-hydrogen recombina- 
tion to form molecules - these events are only allowed in 
the presence of walls. The H2 solid grows layer- by- layer 
from at a rate of 0.5-1 molecular layer per hour. After 
~ lweek, a quartz microbalance revealed a film thickness 
of 150 ± 1 layers^. 



B. Anomalies and experimental results 

Ahokas et alZ observed four anomalies: (1) several 
orders-of-magnitude too fast "Overhauser" relaxation, 
(2) a non-Boltzmann a-b population ratio, (3) satura- 
tion of the a-b spectroscopic line fails to give a 1:1 pop- 
ulation ratio, and (4) recombination rates are extremely 
low. Items 2 and 3 will be our main focus. 

Overhauser relaxation. The c to a relaxation is 
expected to be extremely small at the 4.6T fields of the 
experiments. This can be seen from the small mixing 
angle 9 ps 3 x 10~ 3 , indicated in Fig.Q] The mixing angle 
appears in the states as \a) = cos 8 |4,-f) — sin 6* |t4) , and 
\b) = cos 8 |t 4=) — sin (9 ||t)- Any c-a decay mechanism 
by photon emission is suppressed by 8 a within a > 1. 
Ahokas et alZ observe no such suppression: the c-a line 
decays with a time constant of < 5s, similar to the d-a 
decay time. 

Equilibrium populations. The polarization 

n a - n b 

P = ; (1) 

n a + rib 

characterizes the a and b state population. Assuming 
a Boltzmann distribution n a /rib = exp(A a b/T) where 
A a fc Ri 43mK is the difference in energies between the b 
and a states, one expects p = 0.14 at 150 mK. On the 
contrary, Ahokas et alZ measure p = 0.5. 

When one of the states is depleted the system returns 
to this non-Boltzmann value on a time scale of ~ 50hours. 

Saturation of a-b line. An extremely strong rf field 
should saturate the a:b line, driving the population ratio 
to 1:1, or p = 0. Ahokas et alZ obtained a minimum of 
p = 0.2 at high excitation powers. For sufficiently large 
power this saturated population ratio was independent 
of the applied rf power. 

Low recombination rates. As previously described, 
at 150mK Ahokas et al.'sZ recombination rate is much 
smaller than expected, negligible on a time scale of weeks. 
In contrast, at T = IK their recombination rates are 
consistent with previous studies^. 

Hole burning. Ahokas et al. applied a magnetic 
field gradient and a rf field to saturate the a-b line in a 
millimeter sized region of the sample. The spectral hole 
recovered in a time similar to that in the homogeneous 
case, indicating that the nuclear spin-relaxation is some- 
what faster than spin migration. This would seem to 
indicate that the atomic hydrogen defects are immobile 
on regions much larger than a millimeter. 

III. SCENARIOS 

In this section we evaluate some previously proposed 
scenarios for these phenomena and suggest a new onei. 
For each scenario, we present the idea, examine its con- 
sistency with Ahokas et al.'sZ experiments, consider pos- 
sible microscopic mechanisms, and give testable predic- 
tions. 
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Bose statistics and Bose-Einstein condensation 



Idea 



Ahokas et all suggested Bose-Einstein condensation 
(BEC) as a possible mechanism to explain the departure 
from the Boltzmann distribution. In a BEC the lowest 
energy mode becomes "macroscopically occupied," lead- 
ing to an excess of a-H. Neglecting interactions, the BEC 
transition temperature T c for a homogeneous system of 
spinless particles with density p and effective mass to* is 



G(3/2) 



2/3 



2tt?i 2 

fe^TO* 



(2) 



Here ( is the Riemann zeta function; C(3/2) w 2.61—. 
Including the b states is straightforward and makes only 
small changes: for example if A (, = one would divide 
the density by a factor of 2. Note that as to* becomes 
larger the transition temperature becomes smaller. 



2. Phenomena explainable 

This scenario can in principle explain the non- 
Boltzmann equilibrium ratio rib/n a . It does not pro- 
vide an explanation of the inability to saturate the a- 
b line, the slow recombination, or the fast Overhauser 
cross- relaxation. 



0.2mm hole persisted for T pcrs > 50hours. Our argument 
will relate macroscopic motion (which fills in the spectral 
hole) to microscopic motion (the defect tunneling). We 
assume diffusive motion, where the characteristic time 
between collisions is longer than a tunneling time. 

There arc at least three mechanisms by which the spec- 
tral hole can heal: the excited atoms can spontaneously 
undergo a transition back to the a-state, spin exchange 
collisions can lead to spin diffusion, or a-state atoms can 
diffuse back into that region of space. Neglecting all but 
the last process gives us an upper bound on the atomic 
diffusion constant D. 



D < 



10~ 8 cm 2 /s- 



(3) 



Throughout this argument we will aim to produce an 
order-of-magnitude estimate, and use the symbol "~" to 
indicate that we neglect constants of order unity. This 
diffusion constant can be related to the microscopic col- 
lision time T co u and the mean velocity v by 



D ~ v£. 



(4) 



with I the mean free path. In the effective mass approxi- 
mation, one would expect thermal effects to yield a mean 
velocity v ~ ^/hsT/m*. An lower bound for the mean 
free path is given by the lattice constant I > d ~ 3iU£. 
Thus we arrive at the following lower bound for the ef- 
fective mass 



3. Consistency with experiment 

Transition temperature and densities. As re- 
ported in Ahokas et all, for density p ~ 10 18 cm -3 and 
effective mass to* similar to bare mass m, the ideal Bose 
gas transition temperature in Eq. ([2} is T c ~ 30mK, far 
below the experimental temperature. To address this 
inconsistency Ahokas et all suggest that phase separa- 
tion may concentrate the defects to locally higher densi- 
ties. For example, if the defects phase separated so that 
their density was p~3x 10 19 cm -3 , then with to* = to 
the transition temperature would be T c = 170 mK, 
and one would reproduce the observed ratio of n a /nb at 
T = 150 mK. Such a powerful concentrating mechanism 
would have additional consequences, such as increased 
recombination rates. 

A key question is how large the effective mass is. In 
the following subsection we use the experimental hole- 
burning data to constrain the effective mass, finding that 
it is sufficiently large to completely rule out simple Bose- 
Einstein condensation of defects at the experimental tem- 
peratures. 

Estimate of effective mass. We expect that the 
effective mass to* of the defects is much higher than that 
of the free atoms. Here we bound the effective mass by 
considering Ahokas et aZ.'s measurement of the lifetime 
of a localized spectral hole. They found that a w ~ 



i* > k R T 



10~ 8 cm 2 /s 



10 s toh, 



(5) 



where to# is the bare hydrogen mass and the last ap- 
proximate equality is for the experimental temperature 
T = 150mK. This effective mass is several orders of mag- 
nitude too large to allow BEC at experimentally relevant 
temperature scales, regardless of the defect concentra- 
tion. This argument has neglected interactions and in- 
homogencities: phase separation or some "local" BEC's 
that are uncoupled, thus disallowing global transport, 
would invalidate the arguments leading to our bound. 



4- Microscopic mechanism 

Microscopic estimates of tunneling matrix elements are 
beyond the scope of this work. As previously discussed, 
the exchange reaction mechanism is expected to be the 
dominant pathway^. 



5. Experimental predictions 

Polarization temperature dependence. Perhaps 
the most easily testable prediction of this model is 
the temperature dependence of the polarization, p = 
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(n a /n b - l)/(n a /n b + 1) with 

na Jd 3 kn[e a (k)/(k B T)} 

n b fd*kn[e b (k)/(k B T)} W 

with n(x) = 1/ (e x — 1). The energy dispersion of atoms 
is e a (k) pa h 2 k 2 /2m*, and Q,(fc) « H 2 k 2 /2m*+A ab , where 
A ab is the b-a energy difference. Above the BEC transi- 
tion temperature T c , the integrals yield 

na 93/2 (^ /T ) 

n b ff 3 /2 ( e ^-A Q ,)/T) ^ ) 

where g a {x) = ■ x J / j a is the polylog function and is 
sclf-consistcntly determined to set N, for a homogeneous, 
three-dimensional gas. The same expression holds in di- 
mension d with 3/2 replaced by d/2. Below T = T c , one 
instead finds 

n 6 V r / 53/2 (e A « b/T ) 

While the effective mass sets the density at T c , it does not 
appear in this expression. The polarization depends only 
on A ab /T and T/T c . To produce the observed n a /n b = 3 
at T = 150mK, one needs T c = 170mK. 

Fig. [5] shows p(T) for Boltzmann and Bose condensed 
(assuming T c = 200mK) gases. Accurately measur- 
ing p(T) would clearly distinguish Bose and Boltzmann 
statistics. 

Equation J7]) shows that Bose statistics can affect the 
ratio n a /n b even if the system is non-condensed. How- 
ever, Eq. Q bounds n a /n b < C(3/2)/ ff3/2 (e- A '"'/ T ) = 
2.3, which is insufficient to explain the experimentally 
observed n a /n b = 3. 




Temperature (mK) 



FIG. 2: Fhe polarization versus temperature for the Boltz- 
mann case (solid line) and the Bose-condensed case (dashed 
line), from Eq.0 

Transport. A second signature of BEC is superflow. 
For example, the sample could be incorporated into a tor- 
sional oscillator, providing a measurement of a possible 
nonclassical moment of inertia /: below T c , the superflow 
decouples from the cell, and I decreases. Mounting the 



sophisticated hydrogen growth and measurement equip- 
ment in an oscillator would be challenging, as would the 
difficulty of working with such small samples. 

Bimodal Cold Collision Shifts. Bose Einstein con- 
densation also has implications for the cold collision shifts 
in the atomic spectra. Insofar as the interaction may be 
described by the s-wave scattering length the a-d line will 
be shifted byi£-i£ 

Airh 2 

Suj = g 2 (0) {a n - a u ) (n) (9) 

where g 2 (^) = (ip'' (r)ip^ (0)ip(0)ip(r)} / (n) 2 , (n) is the 
average density, and a-fi, an are the a-a and b-b scat- 
tering lengths, respectively. For a noninteracting BEC 
32(0) = 1 while for a normal gas 32(0) = 2. For a 
partially condensed gas one in fact sees a bimodal spec- 
trum with two peaks: one from the condensed atoms 
and one from the noncondensed atoms. This technique 
revealed BEC in magnetically trapped spin-polarized hy- 
drogen gasi£. 

Thermodynamics and collective excitations. 

The BEC phase transition can in principle be directly 
observed by monitoring thermodynamic quantities such 
as specific heat. Due to the small number of H atoms, 
the signal should be quite small. Similarly, the presence 
of a supcrfluid component would lead one to expect a 
second-sound mode, which could be excited (for exam- 
ple) via localized heating of the sample. 

B. Local Bose-Einstein condensation 

1. Idea 

Next we pursue the idea of "local BEC", where the 
defects aggregate in small disconnected regions, each of 
which contains a condensate, but which have no relative 
phase coherence. 

2. Phenomena explainable 

This model can explain the non-Boltzmann ratio 
n b /n a , and the slow transport observed in the hole burn- 
ing experiments. It does not provide an explanation of 
the slow recombination, failure to saturate, or fast Over- 
hauser relaxation. 



3. Consistency with experiment 

The arguments from Section UlI Al about the polariza- 
tion go through without change. The slow recovery in 
the hole burning experiment is readily explained if the 
disconnected condensates are smaller than 0.1mm. Fur- 
thermore, local clusters are naturally expected if there is 
the dramatic sort of concentrating mechanism described 
in Section InTAl 
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4- Microscopic mechanism 

An attractive interaction between defects can lead to 
clustering. A long distance phonon mediated attraction 
is expected for this system-^. Inliomogeneities in the 
molecular sample, or its environment could also lead to 
clustering. For example, it has been observed that the 
ortho and para molecules phase separate. Furthermore, 
the surface that the sample sits on creates inhomogeneous 
strains, which couple to defects. This can possibly cause 
them to accumulate. 



5. Experimental predictions 

The temperature dependence calculated for global 
BEC is unchanged for local BEC, and one again expects 
a double-peaked ESR spectrum. Jointly observing these 
would provide a "smoking gun" for local BEC. We note 
that at sufficiently cold temperatures the puddles phase 
lock giving a global BEC. However, the transition tem- 
perature would be effectively zero since atoms would have 
to tunnel macroscopic distances between concentrated re- 
gions. The local concentration of Hf , regardless of BEC, 
may be diagnosed by examining the dipolar shift of spec- 
tral lines due to H-H interactions. 



C. Nuclear spin dependent Density-of-states 

Here we outline a non-BEC scenario for the anomalous 
observations. 



1. Idea 

If the degeneracy of the a and b states were g a and gi,, 
one would expect that n a /nb = (g a /gb)e^ Aab ■ Thus if a 
mechanism could be found to enhance (g a /gb), then one 
could explain the observed ratio of n a /rib. Assuming such 
a relative enhancement of the density of states, a strong 
RF field would lead to a saturated ratio (rt a /"-b)sat = 
9a/gb- 



2. Phenomena explainable 

Both the equilibrium n a /nb and the RF saturated 
(n a /nb) sa t can be explained by this model. We also pro- 
vide a scenario whereby the Overhauser relaxation is en- 
hanced. Within this model, the low recombination rates 
could be a consequence of the defects being immobile. 



3. Consistency with experiment 

At strong excitation powers the polarization p = 
("tot, a - nb)/(ntot,a + n b ) saturates to 

Psat = g —^r (10) 

5+1 

where g = g a /gb- Meanwhile, the thermal polarization is 

_ ,gexp(A ab /(fc B r))-l g - 0.75 
Pthcm ~ ff exp(A ab /(fc B T)) + l~.g + 0.75 (U) 

where the last equation holds for Ahokas et aL'si experi- 
ments, where A a f, = 43mK and T = 150mK. If one takes 
g = 2 one finds p sat = 0.33 and pthcrm = 0.45. Ahokas 
et al£ experimentally find p sa t — 0.2 and Ptherm = 0.5. 
One should contrast this with the naive expectation of 

p sat = 0.5 and Ptherm = 0.14. 



4- Microscopic mechanism 

There are very few mechanisms whereby the molecu- 
lar hydrogen matrix can change the degeneracies of the 
atomic hydrogen hyperfinc states. The most plausible 
source would be to consider nuclear spin dependent in- 
teractions with o-H 2 . Such interactions can be produced 
through spin-orbit coupling in the presence of a bias mag- 
netic field. Hybridization of the molecular and atomic 
levels could in principle lead to sufficiently drastic re- 
arrangements of the hyperfinc states to affect their de- 
generacy. Such a strong interaction would presumably 
have other spectroscopic implications, such as a severe 
renormalization of the a-b splitting. Unfortunately, the 
experiments observe that the a-b splitting is changed by 
only 0.1% compared to its vacuum value. 

If there is significant hybridization of the atomic and 
molecular states, then the symmetry which forbids the 
a-c transition would generically be broken. This would 
be a source of the fast Overhauser relaxation. 



5. Experimental predictions 

One would expect that at low temperatures, g should 
be roughly independent of temperature. Thus the tem- 
perature dependence of the polarization in Eq. (|lll) can 
be compared with the experiment. For any given mech- 
anism additional predictions are possible. As one ex- 
ample, in a spin-orbit mechanism hinted at above, the 
physics should also manifest in the populations of the 
o-H 2 states, which can be probed spectroscopically. Fi- 
nally, under further assumptions, one may be able to pre- 
dict the power dependence of the polarization saturation 
experiments. 
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IV. OTHER OBSERVATIONS 

We would like to point out three other possibly rel- 
evant observations. First, Ceperley et al. have shown 
that the surface of small p-H 2 clusters in a vacuum are 
supcrnuid2£. Similarly, Cazorla et al. have shown that on 
small length scales 2D p-H 2 has superfluid correlations^. 
Analogous effects are predicted for defects in solid 4 Hc, 
including grain boundaries, dislocations, and amorphous 
regions^— . Ahokas et al.'s crystals are of very high 
quality, and it is doubtful that grain boundaries and dis- 
locations are playing an important role. On the other 
hand, one could imagine that the the interface of o- and 
P-H2 clusters could play a similar role. 

Second, it seems possible that the spectroscopic 
anomalies are related in some way to torsional oscilla- 
tor mass decoupling observed in solid hydrogen^ at sim- 
ilar temperatures. It is important to note that Clark et 
al. 25 ruled out global supersolidity in their experiments 
by comparing the response of the system in an open and 
blocked annulus. Local supersolidity, however, can also 
give rise to a small period drop, which would be present 
in both the open and blocked annuli. 

A third observation is that magnetic ordering transi- 
tions — for example ferromagnetism — would alter the 
ratio of a-state to 6-state population. One can eliminate 
ferromagnetic ordering, as it would give rise to an observ- 
able energy shift. On the other hand, perhaps a different 
form of magnetic order is playing a role. 

V. SUMMARY, CONCLUSIONS, AND 
CONSEQUENCES 

We reviewed the unexplained phenomena seen in ex- 
periments of Ahokas et al£. Wc enumerated a number of 
possible mechanisms which could be involved in produc- 



ing the observed phenomena. In particular, we gave de- 
tailed consideration to the idea, first introduced in Rcfj£ 
that the non-Boltzmann ratio n a /rib may be due to Bosc- 
Einstcin condensation of atomic hydrogen. We conclude 
that global Bose-Einstein condensation is not consistent 
with other experimental observations. 

Although wc present several other scenarios, we find 
that none of them are wholly satisfactory. Although some 
of the phenomena can be explained by local BEC, it fails 
to provide a mechanism for the unexpected saturation 
population (n a /nb) S at when a strong RF field is applied. 
We do find that all of the phenomena would be consistent 
with a nuclear spin dependent density of states. However, 
we are unable to provide a microscopic mechanism for 
this density of states. 

Ultimately, substantial experimental work will be nec- 
essary to clarify the situation. Our arguments make a 
strong case that measuring the polarization's tempera- 
ture dependence is a promising first step, and suggests 
other experimental signatures — especially in transport 
and spectral features — that would clarify the phenom- 
ena. 

During the preparation of this paper, new results came 
out from Ahokas et al£&, which introduced new myster- 
ies. In particular they observe substantial density and 
substrate dependence of the population ratio n a /nt. All 
of our considerations remain valid, with the additional 
clue that whatever the underlying mechanism is, it must 
involve the surface of the sample, and be sensitive to den- 
sity. For example, the formation of superfluid domains 
could be influenced by the substrate, or magnetic impuri- 
ties in the substrate could interact with atomic hydrogen, 
leading in some way to the unexpected density of states. 
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